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Abstract A time-dependent double-barrier option is a derivative security that delivers the terminal value 
(1){St) at expiry T if neither of the continuous time-dependent barriers b± : [0,T] M.+ have been hit 
during the time interval [Q,T]. Using a probabilistic approach we obtain a decomposition of the barrier 
option price into the corresponding European option price minus the barrier premium for a wide class of 
payoff functions 0, barrier functions b± and linear diffusions (5(),£[o j-j. We show that the barrier premium 
can be expressed as a sum of integrals along the barriers b± of the option's deltas A± : [0, T] ^ M at the 
barriers and that the pair of functions solves a system of Volterra integral equations of the first 

kind. We find a semi-analytic solution for this system in the case of constant double barriers and briefly 
discus a numerical algorithm for the time-dependent case. 

Keywords Time-dependent single- and double-barrier options, local time on curves, Volterra integral 
equation of the first kind, delta at the barrier 
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1 Introduction 

Barrier options play an important role in modern financial markets. They are a less expensive alternative to 
European options and trade in large volumes particularly in foreign exchange. A knock-out double-barrier 
contract is nullified if either of the two barriers is breached by the underlying asset price process during the 
life of the option, and delivers ^{Sj), for some predefined payoff function , otherwise . A knock-in option 
becomes a European option with payoff (j) if one of the barriers is hit by the asset price process before 
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time T, and expires worthless otherwise. Since simple parity relations exist for the prices of knock-in and 
knock-out contracts, we shall concentrate only on examining the later. 

The main result of this paper is given by the following representation formula 

y(0,5o) = (p(0,5o) - ^-e-'-^ y\4t)q,{So,b4t))dt^ l\+it)q,{So,b+{t))dt^ , (1.1) 

where V{Q,So) is the current time-dependent barrier option price, (p{Q,So) is the current price of the corre- 
sponding European payoff and the function is closely related to the transition density of the process St 
(see formula (12.71 )). which is a local volatility process given by the stochastic differential equation (SDE) 
in ( 12.11 ). The functions A^,A+ : [0,T] R can be interpreted as the limiting values of the option's deltas 
^(f,5() as the asset price process St approaches either of the two barriers b-(t),b+{t) at time t (see The- 
orem l2.5l for more details). Note that, since the option's payoff is non-negative, the delta at the lower (resp. 
upper) barrier is positive (resp. negative) making the barrier premium in the above formula negative, as one 
would expect. By Theorems 12.51 and 12.71 the pair of functions exists and solves a system of two 

Volterra integral equations of the first kind, given by ( 12. lit . 

An important feature of representation dl.lb is that it can be used for hedging time-dependent barrier 
options. Once the system of Volterra integral equations in ( 12.1 It is solved (numerically or otherwise), the 
barrier option price can be obtained by computing the one-dimensional integral in formula dl.ll ), which, 
from a numerical point of view, can be done very efficiently. Therefore an entire "spot-ladder" of option 
prices (i.e. a vector of values y(0,5'o), where Sq ranges over a discrete subset in some interval) can be 
obtained with little numerical effort, for we are only solving the system of Volterra integral equations once. 
Moreover, since the function qt is available in semi-analytic form in most models used in practice, spot- 
ladders of deltas and gammas (i.e. vectors with coordinates ^(0,5o), ^^{0,Sii) respectively, where Sq 
takes values in a "discrete" interval) can be found by differentiating formula dl.lb . once we have obtained 
the solution of the system in d2.1 lb . This feature of our pricing algorithm is critical for the risk 

management of barrier option portfolios because spot-ladders are one of the most important tools used by 
traders for understanding their exposure to adverse movements in the underlying market. 

Hedging a down-and-out call, when the barrier level is below the strike, is not dissimilar to hedging the 
corresponding European option, because the presence of the barrier does not destroy the convexity of the 
payoff function and hence the delta ^[t.St) remains bounded throughout the life of the option. In the case 
of a double-barrier knock-out call option the situation is radically different since the barrier option price is 
non-convex close to the upper barrier at any time t before expiry. As mentioned earlier, the value 4+(f), 
where 4+ is the function in (11.11 ). is a good approximation for the delta ^{t^St), when St is close to the 
upper barrier, and can hence be used for hedging. We should stress here that in practice the pair of functions 
(4+, 4-) arises as a solution of a non-singular system of linear equations (see Subsection 13.31 and BTI . 
page 179, equation (4.1)) obtained by discretizing the Volterra integral equations from d2.11l ). rather than 
from a numerical scheme for partial differential equations, where the first derivative is approximated by a 
difference quotient that can be unstable close to the boundary since the value function in that region changes 
at a very rapid pace. 
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Time-dependent barriers arise naturally in financial markets even if the barriers in the option's contract 
are constant. Let St denote the foreign exchange rate and let the functions Rd,Rf ■ [0, T] R describe the 
deterministic term structures of interest rates in domestic and foreign markets respectively. The assumption 
that term structures of interest rates are deterministic is very common in the foreign exchange markets 
as the majority of barrier option contracts are short dated (with maturities up to one year) and have little 
dependence on the stochasticity of interest rates. The forward process :— 5'rexp(— f (/?^(f) —Rf{t))), 
which must be a martingale since it is proportional to an asset divided by the domestic bond, is often 
modelled directly instead of the FX rate St- It is clear that the original barrier option's contract with constant 
barriers translates into a contract with time-dependent barriers for Fi . Furthermore, by modelling the forward 
process directly we can extend the representation in ( 11.11 ) to the case of time-dependent interest rates. Note 
that the process St studied in Section |2] (see ( 12.11 )) has a constant drift. But if the functions Rd,Rf are in 
C^([0, T]), as they usually are since market participants do not like to see kinks in their term structures, we 
can model the forward Ft by dFt — Fta{Ft)dWt and price the equivalent derivative with barriers b^{t) := 
b±{t)exp{-t{Rd{t)-Rf{t))) using (O- 

A feature frequently encountered in barrier option markets is the existence of discontinuous barriers. 
The barriers are usually step functions or simply stop being active at a certain time before expiry. Since 
formula ( II. Il l works for discontinuous payoffs, it can be be applied to non-continuous barriers by "backward 
integration". Time steps would in this case be determined by the intervals of continuity of the barriers. The 
procedure starts at the end of the last such interval where the payoff is known and uses Theorem 12.71 to 
determine the payoff function at the beginning of that interval. This produces an equivalent problem with a 
smaller number of time intervals, so the same procedure can be reapplied until we obtain the option value 
at the current time. 

The key idea in behind the proof of Theorem l2.5l which yields representation dl.ll ). is in some sense 
analogous to that used for finding the integral equation for the optimal exercise boundary in the American 
put problem (see Theorem 4.1 in |31| for a survey account and II2TI for one of the original derivations). 
The smooth-fit principle in the American option problem implies that the value of the first derivative of the 
option price at the exercise boundary is known, which allows us to obtain a non-linear integral equation for 
the exercise boundary by applying Ito's lemma and taking expectations. In the case of barrier options, the 
boundary of the region is specified in advance but the first derivative (i.e. the option's delta) at the barrier 
is clearly unknown. By judiciously applying Peskir's change-of-variable formula (see appendix|A]and ||35]| 
for more details) and taking expectations as in the previous case it is possible to obtain a Volterra integral 
equation of the first kind for the first derivative. Unlike with the American put option, where the integral 
equation is non-linear in the exercise boundary, in the case of barrier options we obtain a linear equation 
which, when discretized, yields an upper-triangular linear system for the unknown function that can be 
solved directly (see Section[3]and 1411 for more details). 

The literature on continuously monitored barrier options is vast and varied. It appears that two general 
approaches have been formed. In the first one, which mainly deals with constant barriers, one tries to find 
a path-wise (i.e. robust) hedging strategy with European derivatives that either uniquely determines or pro- 
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vides an admissible range for the barrier option price. A model-independent approach of this kind is exem- 
plified in Is). In the case of the Black-Scholes model, and more generally for models with symmetric smiles, 
this approach has been applied to a number of path-dependent derivatives including constant double-barrier 
options (see |4| and [|5j). The second approach consists of calculating directly the expectation in a risk- 
neutral measure of the path-dependent barrier payoff. A probabilistic approach using Laplace transforms 
for constant double-barrier call and put options in the Black-Scholes model is described in ifTTl . A method 
using the joint density of the stock, its maximum and its minimum to find the price of time-dependent barrier 
options in the Black-Scholes model was pioneered in |27|. Boundary crossing probabilities for Brownian 
motion have been used in ll32l to price single-barrier options where the underlying asset price process 
has deterministic time-dependent drift and volatility. In Bol it is shown that the time-dependent double- 
barrier option problem for geometric Brownian motion can be reduced to the constant barriers case by first 
transforming the state-space and then time. A static hedge using calls and puts for a time-dependent single- 
barrier option is described in IT]. The result applies to linear diffusions with compound Poisson jumps, but 
the hedging strategy depends on knowing the values of the barrier contract one is trying to hedge at certain 
times before expiry. This deficiency was also noted in [22] (see page 106), where a simplified derivation of 
the main result from HI is given in the case of diffusion processes. More recent work on time-dependent 
double-barrier options for the same kind of asset price process using analytic tools such as Fourier trans- 
forms. Green's functions and complex integration can be found in 8181 . Ill II . and |[T9l and 1341 respectively. 
Spectral methods are applied to find constant double-barrier option prices in the class of CEV models in Js] . 
Laplace transforms and Wiener-Hopf factorization are used in ||24| to obtain prices and Greeks for constant 
barrier options where the logarithm of the underlying asset price process is a generalized hyper-exponential 
Levy process. This class of processes contains VG, NIG, CGMY and other models that are of relevance in 
finance. Chapter 12 in ||29| contains a wealth of analytic methods for pricing a variety of barrier options 
(time-dependent double barriers with and without rebate) in specific modelling frameworks (GBM, CEV, 
Heston) using the theory of partial differential equations. A local time approach has been pursued for the 
study of the static superhedging of barrier options (see L26J ) and the decomposition of European options 
with convex payoff functions (see i6j). 

In this paper we address the question of pricing time-dependent single- and double-barrier options 
where the underlying asset price process is a linear diffusion with mild regularity conditions on its volatility 
function. Our approach is entirely probabilistic, and combines the two approaches discussed in the previous 
paragraph. We do not make use of first passage time distribution, which is prohibitively comphcated in 
the class of models we are considering. Instead we employ a path-wise analysis of the option price, which 
yields representation ( II. 11 1. The paper is organised as follows. Section|2]contains statements and proofs of 
our main results. In Section[3]we propose a semi-analytic solution (using Laplace transforms) of the system 
of Volterra integral equations that arises in Theorems l2.5l andl 2.7l in the case of constant double barriers. We 
also discuss discretization methods for the general time-dependent barrier case. Section|4]considers briefly 
some open questions related to our results and concludes the paper. 
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2 Integral equations for time-dependent barrier options 

In this section our goal is to find the integral equations that characterize the deltas at the barriers and 
consequently the price of any time-dependent barrier option. Before defining precisely the class of exotic 
options we shall consider, let us specify the underlying model that provides uncertainty in our economy. The 
dynamics of the underlying risky security are given by a possibly weak solution (in the sense of Definition 
5.5.1 in ll25l ) of the one-dimensional stochastic differential equation (SDE) 



where the function a : M+ M+ satisfies a(x) > for all x e (0,°°) and is locally Lipschitz continuous 
in the interval (0,°°) (i.e. for any compact set C C (0,°°) there exists a positive constant Kc such that 
\(y{x) — G{y)\ < Kc\x — y\ holds for all x,y e C). These two assumptions are the only regularity conditions 
applied to the function a throughout the paper The constant := r — 5 is the risk-neutral drift given by the 
interest rate r and the dividend yield 5. 

The assumptions on a imply that the volatility function jc i-^ xa{x) is also locally Lipschitz continuous 
in the interval (0, °o) but may vanish at the boundary x = 0. Under these hypotheses Theorem 5.5. 15 in 1251 
yields a filtered probability space {Q,{^t)te[Q.T]^Q)^ with a filtration {•^t)telo.T] that satisfies the usual 
conditions, and processes S — iSt)te[Q.T] W ~ ^t)te^.T\ defined on Q., such that W is a standard one- 
dimensional Brownian motion with respect to {^t)te\^j\ the process S solves SDE ( 12.1b up to an 
explosion time. Furthermore Theorem 5.5. 15 guarantees uniqueness in law of the solution S. 

For some models given by SDE ( 12.1b . the solution S can reach the boundary point zero of the domain 
(0,oo) in finite time with strictly positive probability (e.g. the CEV process, given by a{x)=x^^^ , can reach 
zero if the parameter p is in the interval (0, 1), see [lOj ). In such cases the absorbing boundary condition 
for the process S at zero is assumed in Theorem 5.5.15 of ll25]| . Our aim is to use the measure Q as an 
equivalent local martingale measure for our economy in the sense of Q. The absorbing boundary condition 
at zero is therefore very natural because any other boundary behaviour would in general introduce arbitrage 
(an arbitrage strategy would be to buy the asset when it is worth zero and hold it). 

The solution S of SDE ( 12.1b behaves differently at the other boundary point of its domain. 

Lemma 2.1 The process S does not reach infinity infinite time Q-almost surely. 

For the precise definition of explosion at infinity see 12511 . page 343. Note that Lemma ITT] implies 
that the integrals in ( 12.11 ) are defined on the entire probability space £2 for any fixed time t, because the 
solution process S is Q-almost surely finite during the time interval [Q,t]. For the proof of Lemma IZTI see 
AppendixICl 

A continuous time-dependent barrier b : [0, T] (0,°°) is by definition a continuous function of finite 
variation. In this paper we will mainly be concerned with double-barrier options. In order to define them 
we need two such functions b± : [0, T] (0,°°) which satisfy b- (f ) < b+{t) for all f G [0, T]. For any fixed 
time s £ [0,T] let the stopping time T.s be given by 




(2.1) 



T, :=inf{ve [0,r-i]; S,+^ eR+ - {b^{s + v),b+{s + v)) } 



(2.2) 
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where M+ [0,°°). Note that by definition we have s + T,. < T, {s + <t} e ^, for all t e [0, T], and the 
property s + Ts <t + Tt holds for s <t <T. 

Let / € [0, r] be the current time. By definition the fundamental price Vt of the discounted contract for 
the barrier option with a non-negative measurable payoff function (p : [0,°o) ^ [0, °°) that started at time is 
given by = E [0(5Tfl)/{To=7'}l=^r] where t{TQ=T} is the indicator function on Q (see [23J, Definition 7). 
Since our market is complete, by Theorem 6 in ll23l (see also Theorem 3.3 in [7]), the fundamental price 
of a derivative security is the smallest initial cost of financing a replicating portfolio of that security. 

It was shown in [23] (see Subsection 1.5.4) that the market price of a derivative security equals its 
fundamental price when, in addition to the standard NFLVR assumption of |j9l, we also stipulate the no- 
dominance assumption of Merton [ 301 . No-dominance intuitively says that, all things being equal, market 
participants prefer more to less, and is only violated if there exists an agent who is willing to buy a dominated 
security at a higher price (for the mathematical formulation of the no-dominance assumption see ll23l . 
Assumption 3). No-dominance is shown to imply that there are no bubbles in the price of the underlying 
asset or in the price of a barrier option that is dominated by a call or a put (see 1231 , Proposition 1, and 
Lemma 8 and Theorem 7). The assumption is consistent with a subclass of models given by SDE ( 12. It , 
namely those that have an equivalent martingale measure. For example in the CEV framework ((y{x) = 
xP^^, where p e (0, 1]), which is known to have a unique equivalent martingale measure (see ifTOl ). the 
no-dominance assumption can be made and the fundamental price given by Theorem 12.71 is the market 
price. However the no-dominance assumption cannot be made when the discounted asset price process 
(exp(— /if )5,),£[o.7'] is a strict local martingale (i.e. there is a bubble in the underlying economy), which 
has been shown to be the case for some of the models in our framework (see ll30l and lfT3]| ). In this case 
the market price of the derivative can exceed the fundamental price given by Theorem l2.7l In other words 
the market price of the derivative is strictly larger than the price of the replicating portfolio and little can 
be said about its dynamics (see [|23ll . Subsection L5.3, Example 5). It is not an easy task to obtain general 
necessary and sufficient conditions for the existence of an equivalent martingale measure for the solution of 
SDE ( 12.11 ), a topic that merits further research. In this paper the process V, denotes the fundamental price of 
the barrier option in the economy given by ( 12. It , referred to simply as the "price" in all that follows^ 

Our aim is to find the price V, at any time f e [0, T] of a time-dependent double-barrier contract initiated 
at time zero. In order to do this we consider the process 

Z, :=£[(/) (5,+„)Vt,=nl^']' (2-3) 

which equals the discounted value of an equivalent time-dependent barrier contract initiated at time t . Unlike 
Vt, which is a martingale under the pricing measure Q, the process Z, is not a discounted price process of a 
security in our economy, since at each time t it represents the price of a different security, and hence need 
not be a martingale (see Lemma |272] |. This is somewhat similar to the well-known observation in interest 
rate theory that the short rate (i.e. the rate at which funds can be borrowed for an infinitesimal period of time 
- also known as the instantaneous interest rate) corresponds to a different asset at each time t and therefore 

' Thanks are due to the anonymous referee for raising the issue of bubbles and their implications for the pricing of options. 
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need not satisfy any no-arbitrage drift restrictions. Unlike in the case of the instantaneous interest rate, 
the drift of Zf can be determined uniquely and, as we shall soon see, contains all the information needed 
to obtain the current price of the barrier option. Before exploring some basic properties of the process Z 
in the next lemma, note that definitions i2.2\ and (I2.3l l also apply to single-barrier options with obvious 
modifications. 

Lemma 2.2 (a) Let the times s,t € [0,T] satisfy s < t. Then the inequality W^t\-'^s\ > holds almost 
surely in (i2 , Q). If either the upper harrier b+ is present or the random variable {St ) is in L} [Q. , Q), 
the process Zt is a non-negative submartingale. 

(b) Assume that the payoff function (j) ; M+ ~* R+ is continuous on the complement of a finite set of 
points where it is right-continuous with left limits, and that, if b^ is not present, the payoff (j) (St) 
is in L^{Q,Q). Let the log-normal volatility O be locally Lipschitz continuous in the interval (0,oo) 
and assume that it satisfies o{S) > Q for all S G (0,°°). Then the process Zi has a continuous mod- 
ification of the form Zt — Z{t,St), where the continuous function Z : [0,T] x IR+ M+ is given by 
Z{t,S) %,s [<p{S,+,,)I{,+,,=T}\- LetC:^ {(f,5) G [0,T) x M+; <S< b+it)}, B+ {{t,S) € 
[0,r) X M+; S > b+{t)} and B := {{t,S) G [0,r) x M+; S < b-{t)} be open subsets of the domain 
[0,T) X R+. Then Z vanishes on the set B UB+, is of order C''^(C) and satisfies the partial differ- 
ential equation Zt{t,S) + liSZs{t,S) + ^-^-^Z5s(f ,5) = for all {t,S) G C with terminal condition 
Z{T,S) = <j){S) for S Cz {b-(T),b^{T)) and boundary conditions Z{t ,b±{t)) = for all t Cz [0,r]. The 
same is true for a single-barrier option price with appropriately modified boundary conditions. 

From now on we shall assume that we are working with the modification of the process Z, given in (b) 
of Lemma 122] i.e. we will assume that the process Z, is a continuous submartingale. Note also that the 
statement in (a) of Lemma 12. 2l is intuitively clear If the underlying asset price is between the barriers, the 
process Zt is a true martingale up to the first time St hits a barrier, because before that stopping time Zt 
equals the discounted barrier option price. Since Z, is non-negative, if it were a martingale it would have to 
stay at zero from that moment onwards. But as soon as the stock returns to the interval between the barriers, 
the process Z, assumes again a strictly positive value. Such behaviour makes its mean drift upwards with 
time. We will now give a straightforward but rigorous proof of this fact. 

Proof. Pick s,t G [0,T] such that ^ < f . Note that i + T.^ < f + T, for all paths in £2 and therefore we have the 
inclusion {s + = T} C {f + Tf = T} and the identity I{s+t:s=t} — t{s+Zs>t}t{t+z,=T}- We can now rewrite 
Zs, using the tower property and the fact {i + > f} G in the following way 

Z, = Em{St+r,)I{s+r,>,}I{r+r,=T}\M^^] = E[Z,/{,,>,_,} |^,] < E[Z,|^J. 

The last inequality holds because (j), and hence Z,, is non-negative. If either of the two integrability con- 
ditions in (a) of Lemma |Z2l are satisfied, we get E[Z(] < °o for all t G [0,r], which implies that Z, is a 
submartingale. This proves (a). Part (b) in the lemma is a well-known fact about barrier options. It suffices 
to note that the statement (b) is a special case of Theorem lB.ll in Appendix IB] r-i 
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Part (b) of Lemma l272l implies that the partial derivative Zs is a continuous function in the open region 
C, but the lemma says nothing about the behaviour of Zs at the boundary of C. A key step in obtaining 
the integral representation for the double-barrier option price (see equation ( 12.101 ) in Theorem 12.51 ) will 
be the application of Theorem lA.il to the function Z : [0, T] x R+ IR+ given in (b) of Lemma l272l This 
step requires a certain regularity of the function Z and its first derivative Zs close to the boundary of C. 
In principle the limit of the delta of the double-barrier option price may not exist as the underlying asset 
St approaches the boundary of the region C. It does not come as a surprise that additional hypotheses on 
the regularity of the payoff function ^ : R+ M+ as well as of the barriers b± : [0, T] M+ are required 
for the function Z to satisfy the assumptions of Theorem lA.il Lemma [23] gives sufficient conditions for 
functions (j) and b± that guarantee that the first spatial derivative of the solution Z of the PDE problem in 
(b) of Lemma 122] does not blow up at the boundary of the region C. 

Lemma 2.3 Let the continuous barriers b±: [0,T]-^ M+ be twice-dijferentiable and assume that the payoff 
function (p : [b-{T),b+{T)] M+ satisfies (p{b- (T)) — (p(b+{T)) — and is twice-differentiable with the 
second derivative (j)" : {b-{T),b+{T)) — > R which is Holder continuous of order a € (0, 1). If only the 
lower barrier b- is present, we additionally assume that the random variable ^{St) is in L^{Q,Q). Then 
the following holds. 

(a) The limits 

A+{t):^\imZs{t.b+{t)-e) and A^{t) ■.^\imZs{t,b-{t) + e) 

£\0 £\0 

exist for all t e [0, T] and are uniform on the interval [0, T]. 

(b) For some 5 > we have supQ^^^gV {Z{- ,b+{-) — e)){T) <oo and supQ^^^gV {Z{- ,b^{-) + e)){T) < °°, 
where V{g){T) denotes the total variation of a function g ; [0, T] — )■ M. 

Properties (a) and (b) hold in the time-dependent single-barrier case with obvious modifications. 

Lemma [23] is a consequence of Schauder's boundary estimates for parabolic partial differential equa- 
tions. For the proof see Appendix |C] Note also that the uniform convergence in the lemma, together with 
Theorem lB.il implies that the delta at the barrier A±{t) is a continuous function of time for all t e [0, T], if 
the barrier functions b± and the payoff satisfy the assumptions in Lemma lZJ] This should be contrasted 
with the known behaviour of the delta of an up-and-out call option which goes to minus infinity if, close to 
expiry, the underlying asset approaches the barrier level. 

The task now is to understand the path- wise behaviour of the process {Zt)te[(i.T] - For this we will need the 
important concept of local time. Recall that the local time (at level a G R) of any continuous semimartingale 
X = (^r)(e[o.7'] on the probability space (X2,Q) can be defined as a limit 

U;{X) := lim i I" I[,,,+,){X„)d{X,X)u 
£\o e Jo 

almost surely in Q (see f3S\, page 227, Corollary 1.9), where {X,X)t is the quadratic variation process as 
defined in Ii38j . Chapter IV, Theorem 1.3 and Proposition 1.18. Notice that this definition can be easily 
extended to a local time of X along any continuous curve b: [0, T] — > R with finite variation by Lf{X) := 
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Lf (X — b), since the process X — b is still a continuous semknartingale and the equality {X — b,X — b)t = 
{X,X), holds for all t. For times < f < v < T we denote the local time of X between t and v by L"^,{X) : = 
L"{X) — L"{X). It is well-known that the map t ^ L"{X) is almost surely a non-decreasing continuous 
function. With this non-decreasing process one can associate a random measure dLf{X) on the interval 
[0,r] the support of which is contained in the set {t E [0, T]; X, = a} (see |38J, page 222, Proposition 1.3) 
By Lemma |2!2] we know that the process Z = (Z,),g[o j-j is a non-negative continuous semimartingale. 
Therefore we are at liberty to apply the Tanaka formula (see ll38l . page 222, Theorem 1.2) at level to the 
non-negative process (Z(),£[o.7']' thus obtaining the following path- wise representation 

Z,=Z,+ J^' I^z^^^o}dZ, + Ul,{Z) for 0<f<y<r. (2.4) 

Using the representation in (I2.4l i we can prove the following proposition, which will play a central role in 
all that follows. 

Proposition 2.4 Assume that the payojf function (j) and the barriers b± satisfy the assumptions ofLemma \2.3\ 
and let t,v be two elements in the interval [0, T] such that t <v. If the upper barrier b+ is present, the pro- 
cess {Jj I{z„>Q}dZii)y^^t.T] i'> ^ continuous martingale and hence representation \2.4^ is the Doob-Meyer 
decomposition of the submartingale (Zy)^,£[; j-j. The following equality must therefore hold almost surely 

Z, = E[Z„|^,] - iE[LO,(Z)|^,]. (2.5) 

Assuming that ^{St) is in L^(X2,Q), the representation i2.5i holds also in the case where only the lower 
barrier b- is present. 

If time V equals expiry T and time f equals the current time, the equality in Proposition 12.41 yields a 
representation of the barrier option price at f as a sum of the current value of the European payoff Zj 
and the expectation of the local time from now until expiry. The former quantity is usually available in 
most models in a semi-analytic closed form and the latter will be obtained in Theorem |2.5| bv applying the 
change-of-variable formula from lf35l . Note also that intuitively the stochastic integral J^' 7{2„>o}'^Z„ is a 
martingale because the inegrator Z„ equals a discounted double-barrier option price on the set {Z„ > 0}, 
which is a martingale. 

Proof. Let C denote the domain between the barriers as defined in (b) of Lemma 12.21 Recall that Z,, = 
Z{v,Sy) where the function Z .C ^ M+ is the solution of the PDE in (b) of Lemma l2.2l By Lemma lZD we 
are at liberty to apply Theorem lA. 1 I to the function Z. In differential form we obtain 

dZu = I[h_(^i,)<s,<b+{u)}Zs{u,Su)Su(y{Su)dWu + ^ (jis„=b^iu)}^-{u)dL''- {S) ~ Is^s„=b+{u)}^+{u)dL'l+ (5)) , 

where Z^ denotes the first derivative of Z with respect to S. The inclusion {Zj, > 0} C {b- (m) < 5„ < fe+ (m) }, 
for all u E [Q, T], follows from definition (12.31 1 and therefore implies 

^ t{z„>o}dZu ^ I{z„>o}Zsiu,Su)Sua{Su)dWu. (2.6) 
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The function Z5 is bounded on the domain C by Lemma 12.31 which impUes that the stochastic integral 
on the right-hand side is a continuous martingale starting at zero. By taking expectation on both sides of 
equality (12.4b we conclude the proof in the double-barrier case. However, the same argument can be applied 
if only the upper barrier b+ is present. This is because the integrand on the right-hand side of ( I2.6l l is still 
bounded, making the stochastic integral in (12.6b a true martingale. 

In the single-barrier case with only present, the argument above does not work because the integrand 
in ( 12.61 ) is no longer necessarily bounded. However the same reasoning shows that the identity in ( 12.5b holds 
for the stopped process Z,.at„ = Z(v A T„ , 5vat„ ), where the stopping time T„ is the first passage time of the 
diffusion S into the interval [n, °o) after time f , because the integrand in ( 12.6b is bounded. Jensen's inequahty 
for conditional expectations and the definition of the process Z given in ( 12.31) imply the inequality 

max{E[z2|^,],E[z2^,J^,]} < E[(^)(5r)2|^r]. 

Our assumption on ^ [St ) implies that Z,, and ZvAt„ are elements of the space (i2 , Q) for all large natural 
numbers n. 

By Lemma im we have that lim„^ooT„ is infinite almost surely in X2. In other words we have almost 
sure path-wise convergence lim„^ooZ,,AT„ = Z,,. The Cauchy-Schwartz inequality implies the following 

E[|Z, - Z„A.„ I l^r] - E[/,„<,|Z, - Z„AT„ I m 

< IE[/.„<v|^r]'/'E[|Z„ -Z,.A.„ I V^]'/' 
<2E[/,„<„|^,]i/2E[0(5r)Vf]'/'- 

Since the sequence E[/t„<,.|^(] converges to zero Q-almost surely as n goes to infinity, we obtain 

E[Z,.|^,] = lim E[Z,.A.„ \n =Zt-\ lim E[L«,.^,,_ (Z)|^,] = Z, - iE[L«„(Z)l^r], 

where the last equality follows by the monotone convergence theorem. This concludes the proof of the 
proposition. □ 

Before we proceed to our main theorem recall that, for any point x e {Q,°°) and time t G {0,T], the 
density p{t;x, •) : {Q,°°) — > M+ of the transition function of the underlying asset price process S, given by 
the SDE in ( 12.1b . is characterised by the identity Qx{St & A) = p{t;x,y)dy, where A is any measurable 
set in (0,00). The function p{t;x, •) : (0,°°) R+ is non-negative but does not necessarily integrate to one 
because the process can reach zero (and stay there) in finite time. The existence of p{t;x,y) can be deduced 
from II20II . Section 4. 11, where it is shown that the transition function of a diffusion is absolutely continuous 
with respect to the speed measure (see ll20l . page 107, for the definition of the speed measure). In the case 
of the process S given by ( 12.11 ), the speed measure is absolutely continuous with respect to the Lebesgue 
measure on the interval (0,°°) and hence the existence follows. Furthermore it is proved in ||20l (page 149) 
that the function p{-',-,y) '■ (0,°°) x {Q,T] R+ satisfies the parabolic PDE in (b) of Lemma 12.21 for any 
y G (0,00). This fact will play a crucial part in the proof of Theorem l2.5l (cf. proof of Lemma ( 12.61 )). Note 
also that sufficient conditions for the existence of densities of solutions of one-dimensional SDEs, which 
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are jointly smooth in all three variables, are given in Il39l . This stronger result requires a volatility function 
that is uniformly bounded away from zero and is therefore not suited to our purpose. For most models that 
are of relevance in mathematical finance the densities can be obtained either in semi-analytic closed form 
(see for example ( 12.81 ) and ( 12. 9t ) or numerically. 

The kernels of integral operators appearing in Theorems 12 . 5 1 and | 2 .7 1 are related to the transition function 
of the asset price process {St)te[o.T] ™d will now be specified precisely. The quadratic variation {S,S)t is 
a continuous non-decreasing adapted process and as such defines, for each path of St, a measure d{S,S)t 
on the interval [0,r]. Since the asset price process 5, is a solution of the SDE in ( 12.1b . this measure is 
absolutely continuous with respect to the Lebesgue measure on [0, T] and the Radon-Nikodym derivative is 
given by d{S,S)t = Sf<7{St)^dt. The function qt{x,y) that appears in the kernel of the integral operators in 
Theorems 12 . 5 1 and 12 . 7 1 c an be defined as 

.d{S,S), 



q,{x,y) := p{t;x,y)- 



dt 



(2.7) 

St=y 



where p{t\x,y) is the density defined above. In the case of the geometric Brownian motion we have the 
formula 

where the drift equals jj. — r — d and is the constant variance. The function p{t;x, •) : (0,oo) M_|_ in 
the case of GBM is a true probability density function because the process cannot reach zero. In the case 
of the CEV model, given by ( 12. Il l with absorbing boundary condition at zero and the log-normal volatility 
function a{x) = GqxP^^ where p e (0, 1) and Ob G (0,°°), we have the following closed form expression 
for the function qt: 

q,ix,y) = 2(Toy''(l ~p)fc'/P-2P)(Xyi-4'')i/(4-4p)exp(-X-F)/i/(2_2p) {2Vxy) . (2.9) 



This expression is a consequence of ( 12.71 ) and the formula for the transition density pt, which can for 
example be obtained from Theorem 3.5 in ifTOl . The function z Ia{z) is the modified Bessel function of 
the first kind of order a and the parameters in (I2.9l l are given by 

2(T2(i_p)(exp(2fAl(l-p))-l)' 
X :=kx^^^-Phxp{2ti^{l-p)), 

Y:^ky^^'-P\ 

where jj, is the drift in SDE (12.11 ). We now state one of our main theorems. 



Theorem 2.5 Let St be the underlying process given by ( 12.71 ) and let Zt — Z(t,St) be the discounted price 
of a time -dependent single- or double-barrier option contract, starting at the current time t, given in ( I2.3I ). 
Assume further that the barriers b± : [0, T] — s- R+ and the payoff (j) : — > satisfy the assumptions 
of Lemma 12.31 and that the local volatility function x i-^ (y{x), x € M+, satisfies the assumptions in (b) of 
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Lemma \2?2\ In the case where only the lower barrier is present, we assume in addition that the variable 
<^{St) is in L?(Q,Q). Let (p{t,x) := £(,^[^(57')] denote the discounted current price of the European con- 
tract starting at time t, conditional upon the asset price St being at level x, and let the function qt{x,y) be 
as in ( 12. 7I ). Then the following integral representation for the time-dependent double-barrier option price 
holds 

Z(0,5o)-(p(0,5o)-^ r A^{t)q,{S(,,b-{t))dt + ]- r A+{t)qr{So,b+{t))dt, (2.10) 
2 Jo 2 Jo 

where A± (t) is the limiting value of the delta of the double-barrier option price at b± (t) as defined in (a) of 
Lemma IZil Furthermore the continuous functions : [0,r] M satisfy the following linear system 

of two Volterra integral equations of the first kind 

(p(t,b+(t)) , ^ , , , , , 




\(p{t,b-{t)) 

where the matrix Q{t,u), for <t < u <T, is given by 

(-qi>-t{b+{t),b+{u)) qu-,{b+{t),b-{u))\ 
\-qu-,{b-{t),b+{u)) qu-,{b-{t),b-{u)))' 
In a time-dependent up-and-out ( resp. down-and-out) single-barrier case, representation ( 12.701 ) contains a 
single integral along ( resp. b-). The integral equation that determines the function 4+ ( resp. A^} in the 
up-and-out ( resp. down-and-out) case takes the form of the Volterra equation of the first kind with ± equal 
to + (resp. —): 

1 r'^ 

(p{t,b±{t))±-J^ qu-,{b±{t),b±{u))A±{u)du = Q. (2.13) 

Theorem 12 .5 1 yields an integral representation for the double-barrier option price for a wide variety of 
local volatility models, any pair of time-dependent barriers and any payoff function that satisfy the assump- 
tions in Lemma lZJ] Rather surprisingly, knowing the values of the delta at the barriers for all future times, 
as well as the current price of the corresponding European derivative (recall that {T)) = (p{b+{T)) = 
for payoffs satisfying the assumptions in Lemma l23T l, is enough to obtain the current value of the time- 
dependent barrier option. Note also that both integrals in equation (12.101 1 are negative since 4- (f ) > (resp. 
A^{t) < 0), which intuitively follows from the fact that the barrier option price is increasing (resp. decreas- 
ing) as the asset price moves away from (resp. approaches) the lower (resp. upper) barrier As expected this 
makes the barrier option cheaper than its European counterpart. Representation ( |2.10t therefore decomposes 
the double-barrier option price into the European option price and the barrier premium. 

In order to include the payoff functions (j) that are of interest in applications (e.g. the up-and-out call 
option payoff 0(5) = {S — K)^I(^Qi,^(^j-j){S) or the payoff of a double-no-touch 0(5) = I{b_{T).b+{T))iS)), 
we must relax the smoothness requirements for the function stipulated in Lemma |23] This will be done 
in Theorem 12.71 where we will show that the integral representation for the price (12.101 1 and the integral 
equation for functions A± ( 12.1 lb continue to hold. 

Before proceeding to the proof of Theorem l2.5l we need the following lemma that bounds the growth of 
the function q, defined in |2.7l over short time intervals. 
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Lemma 2.6 Let K be a compact interval contained in (0, °°). Then there exists a positive constant Ck such 
that the inequality 

qu-t{x,y) < 



\/u — t 

holds for all t < u <T andx,y S K. 

The proof of Lemma l2.6l is contained in AppendixiBl Note that the constant Ck in Lemma (12.6b depends 
only on the compact set and the inequality therefore holds uniformly on K. If the function a in SDE (12. Il l 
were uniformly bounded away from zero, the estimate in Lemma i2.6i would hold on the entire domain 
(0-). 

Lemma ( 12.61 ) implies that integral equations ( 12.111 ) and ( I2.13l l have weakly singular kernels and that the 
inequalities q,i^t{b±{t),b±{u)) < J^_^ hold for all m G {t,T] where M is a positive constant (± denotes 
either + or — ). The linear operator in (12.13b (resp. ( 12.111 )) is compact on the Banach spaces of continuous 
functions C([0, T]) (resp. C([0, T]) x C([0, T])) with the supremum norm, and as such has in its spectrum. 
Note that by construction equations ( 12.111 ) and ( 12.13b have a continuous solution. The uniqueness of this 
solution is a much more subtle question, equivalent to asking whether in the spectrum of the operator is 
an eigenvalue. Since equations ( 12.111 ) and ( 12.131 ) are of the first kind and the Fredholm alternative (which 
provides a general answer to the question of uniqueness of solutions for the integral equations of the second 
kind) cannot be used, it is difficult to answer the question in general. However for a time-dependent single- 
barrier case in the Black-Scholes model see Proposition |2.8l Let us now proceed to the proof of Theorem |2.5l 

Proof. Let us start by considering a time-dependent double-barrier option. Let C be the domain between 
the barriers as defined in (b) of Lemma l272l We begin by applying Theorems I A. 1 1 and IB . 1 1 to the process 
Z, = Z{t,St), where the function Z is the solution of the PDE from (b) of Lemma l272l For any pair of times 
f , V e [0, r], such that t < v, we therefore obtain the following path-wise representation 

Z, = Z{t,Sr) + j/"' /{,_(„)<5„<i^(„)}Zs(M,5„)5„CJ(5„)«'W„ + ^ A4u)dLt- (S) - ^ J\+{u)dL^,+ (S), 

where the functions A+ and Zi_ are defined in Lemma [23] The random measures dLii^{S) are by defi- 
nition equal to the well-defined random measures dL^X'^ — b±) and the functions A± are continuous by 
Theorem IB. II and (a) of Lemma 123] and are hence Borel measurable. Since the function Zs : C ^ IR+ is 
bounded, this equality yields a Doob-Meyer decomposition of the submartingale (Zv),,g[( j-j. Since such a 
decomposition is unique, Proposition l2.4l implies the following identity for the finite variation processes 

lUz) = \f^A^{u)dLi-{S)-\j^'A+{u)dL>:;{S). (2.14) 



The main idea for the proof of Theorem 12.51 is to use the equality in Proposition 12.41 to obtain the 
representation of the option price and the integral equations in the theorem. We must therefore find the 
expectation E; 5, [LJ'^,(Z)] using identity (12.14b . Let us start by proving the following. 
Claim. For any continuous function / : [0, T] — ^ R of finite variation and for all f, v £ [0, T], such that t <v, 
the equality 



f{u)dLl^{S) 



f{u)qu-t{St,b±(u))du 
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holds, where qii^t{x,y) is given in (I2.7l i. 

Recall that, since the process {S — i'±),G[o,r] ^ continuous semimartingale, there exists a modification 
of the local time Lf ^,{S — b±) such that the map a i-^ L" ^,{S — b±) is right-continuous and has left limits for 
every v S [t,T] almost surely in Q. The function is therefore Lebesgue measurable and the occupation times 
formula (see fi38J . Chapter VI, Corollai-y 1.6) implies //'/p £)(5„ — ^±(M))t/(5,5)„ = Jq Lf^^,{S — b±)da for 
e > 0. By taking expectations and dividing by e on both sides of this equahty we obtain 

f -ErA [l[o.e)i^u~b±{u))Sl(yiSu)^] du^- /"'e,,^, [l;',v(5-^±)] da. (2.15) 
Jt £ £ Jo 

The integrand on the left equals ^ JI^^I"-!^^ qu-t{St,y)dy and in the limit e \ we obtain qu-t{St,b±{u)) 
for all u £ (f,v]. Since for all small e we have the inequality q„^tiSt,y) < -^^^ for some constant M and 
y e [b±{u),b±{u) + e] (see Lemma IZ6] ), we can apply the bounded convergence theorem to the left-hand 
side of (12.151 1 to obtain jj qu-t{Snb±{u))du for all values St (including St = b±{t)). 

The right-hand side of (I2.15l l will converge to E; 5, ^.{S — b±)\ by the fundamental theorem of cal- 
culus, if we can show that the function a 1— > E,^, \L"^,{S~b±)\ is continuous at a = 0. Tanaka's formula 
yields the following representation for local time 

^L,"y{S-b±)^*I',,v{a)+ I{s„<b±(u)+a}Su(y{Su)dWu + Ils„<b±{u)+a}{lJ-Su - b'±{u))du, 

where 'fr.v(fl) ■— {a — {Sy — b±{v)))^ — [a — [St — b±{t)))'^ and, as usual, (x)+ := max{x,0} for any x G M 
(see |[38l . Chapter VI, Theorem 1.2). Note that the variable ^^a) is Lipschitz continuous in a with a 
Lipschitz constant equal to 1 for all elements in 12. By taking expectation on both sides we find 



%,s, [Ll,iS-b±)\ =2E,,s,['f^,,(fl)]+2E,,s, 



hs„<b±(u)+a} {llSu-b'±{u))du 



(2.16) 



since the integrand in the stochastic integral is bounded and hence the martingale term vanishes in expec- 
tation. The quantity E, ['f{.v(fl)] is continuous in a, while the second expectation on the right-hand side 
can be rewritten using Fubini's theorem in the following way £ F{a,u)du, where the function F is given 
by F(a,M) := Cq{u) + f^^^"''^" {pLy — b'^{u))p{u — t;St,y)dy, the function p denotes the density of the asset 
price process S in the interval (0,°°) and Co(m) '■— —b±{u)Qsi{Su = 0) is a function independent of a. The 
estimate p{u — t;x,y) < for a positive constant C independent of x,y in a compact subset of (0,°°) 

(see Lemma IZSt . implies that the function a F{a,u) possesses a partial derivative that is bounded in the 
following way | ^(a,M)| < for all m e [f , v]. D is some positive constant independent of St- Lagrange's 
theorem now implies that the integral j^ ' F{a,u)du is a continuous function of a. We have therefore shown 
that ( 12.16b is continuous in a and hence proved that the key identity 

Er.s, [Lf_J(5)] =^ qu-t{St,b±{u))du 

follows from (|2.15l l upon taking the limit e \ 0. 

For every path (O in the probability space X2, the function v Lf^(5')((») has finite variation and is 
continuous. Since the same is true of the function / in our claim, we can use the integration by parts formula 
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to obtain the equality J^ ' f{u)dLu^ (S) = /(v)Lf *(5) — j, ' L^^{S)d fu, where df„ is the Radon measure on 
the interval [t,v] induced by /. By taking expectations on both sides of this identity and applying Fubini's 
theorem to the integral on the right, which is justified since local time is a non-negative function, we obtain 
the following sequence of equalities: 



f{u)diJ:^{s) 



= f{v)%.s, [Lf,J (5)] - j[ %s, [L^^u {S)]dfu 
= f[v)%,s, [if,? {S)] - f df, J" qs-r {S, , b± (s) )ds 
= /(v)E,s, [Lf,J (5)] -^"(/(v) (5„/.± (5) 
f{s)qs-t{St,b± {s))ds. 



The third equality follows by Fubini's theorem and the last one is a consequence of the formula for the 
expectation of local time. This proves the claim. 

In order to apply the claim to the identity in ( 12.14b . we need to approximate the continuous functions A± 
on the interval [t,v] by sequences of uniformly bounded continuous functions (/^ : [f , v] K)«gn with finite 
total variation (since the functions A± are bounded on [t,v], we can take piecewise linear approximations 
on a uniform grid in [t,v]). For each path co E Q the dominated convergence theorem implies that the 



equality J,'' A±{u)dLt^{S){co) = lim„^» /;7±(M)t/L°±(5)(a)) holds. Since the functions /jf are uniformly 
bounded by some constant K, the random variables | {u)dLt^ (S) \ are bounded by A'if^'^ (S), which is 
an integrable random variable. Another application of the dominated convergence theorem and the above 
claim therefore yield the following equalities 



b±, 



A±iu)dL':,^iS) 



lim E,^s, 



f^{u)dLi^{S) 



A±(u)qu-t{S,,b±{u))du, 



for any pair of times f , v £ [0, T] that satisfy t <v. 

We can now apply the last equality to equation ( |2.14| i to find the expectation of the local time of the 
time-dependent double-barrier option price. In other words, by Proposition 12 .41 we have the following rep- 
resentation for the expectation of the double-barrier option price process 



E,,5,[Z,,]-Z(f,5,) 



A^{u)qu-tiSi,b^{u))du- 



A+{u)qu-t{St,b+{u))du (2.17) 



for all t,v £ [0,T] satisfying t < v and all values of St- The representation of the double-barrier option 
price ( 12. lot in the theorem can be obtained by taking f = and v = T in equation ( |2.17t . The system 
of integral equations ( 12.1 11 1 for {A+,A^) also follows from formula (I2.17l i by taking v — T, St = b+{t) 
and St — b-{t), and observing that Z{t,b^{t)) ~ Z{t,b+{t)) = for all t e [0, T], since the double-barrier 
contract that starts at the barrier is worth zero by definition. This completes the proof of the double -barrier 
case. The single-barrier case can be obtained by making straightforward modifications to the preceding 
proof. Q 
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Our next task is to relax the assumptions on the smoothness of the payoff function <j) : M+ K+ made 
in Theorem l2.5l This is crucial, as we would like to be able to apply our methodology to the payoffs that 
arise in practice, such as the up-and-out call option payoff (S) — {S — K)'^I{o,b^{T)) (S) or the payoff of a 
double-no-touch (j){S) = I(h-(T).h+(T)) {^)- The following theorem allows us to do precisely that. 

Theorem 2.7 Let (j) : M+ — s- be a payoff function that is continuous everywhere except at a finite set of 
points where it is right-continuous and has left limits. Assume further that the barriers b± : [0, T] — > R+ 
satisfy the assumptions ofLemma \2.3\ and that the asset price process {St)te[0,T] given by ([2.1i . In the case 
where only the lower barrier b^ is present, we assume in addition that the variable (j) (St) is in {£2 , Q). Let 
Z, ~Z{t,St) be the discounted price of a time-dependent single- or double-barrier option contract, starting 
at the current time t, given in ( 12.31 ), and let (pit,x) := M,t,x[<j>iST)I(b-{T).b+{T)){ST)] denote the discounted 
price of the European contract at the current time t conditional upon the asset price St being at level x. 
Then there exist measurable functions : [0,r] M, which are in L'([0,r],OT±(t/f)) and are not 

necessarily continuous or bounded, such that the following integral representation for the double-barrier 
option price holds 

Z(0,5o) = (p(0,5o)-i f A^{t)qt{S(,,b-{t))dt + \ r A+{t)qr{So,b+{t))dt. (2.18) 
2 Jo 2 Jo 

The measure m± is absolutely continuous with respect to the Lebesgue measure and the Radon-Nikodym 

derivative is given by — qt{b±{0),b±{t)), where ± is either + or — . Furthermore the functions A+,A- 

satisfy the linear system ofVolterra integral equations of the first kind given by ( 12.771 ). In a time-dependent 

up-and-out (resp. down-and-out) single-barrier case there exists a measurable function A+ (resp. A^), 

which is contained in L^ {[Q,T],m^{dt)) (resp. L^ {[Q,T],m-{dt))) and is not necessarily continuous or 

bounded, such that the discounted option price Z(0,5o) has the following integral representation 

Z(0,5o) = (p(0,5o)±- / A±{t)qriSo,b±{t))dt, 
2 Jo 

where ± equals either + or —. The integral equation satisfied by the function A^ (resp. A^) takes the 
form dZTil ). 

At first glance Theorems l2 . 5 l and [2 . 7l look similar. The difference lies in the fact that Theorem l2.7l applies 
to a much wider class of payoff functions (j) that do not satisfy the hypothesis of Lemma l273] and furthermore 
invalidate its conclusions. This makes it impossible to apply the key local time formula from Theorem lA.il 
which provided the core of the proof of Theorem 12.51 These analytical difficulties will be circumvented 
by a careful approximation argument yielding the existence of functions A^,A+, which satisfy integral 
equation ( 12.111 1 and give the desired representation for the time-dependent barrier option price. 

As an illustration of the difference between Theorems l2.7l and [2.5l consider the following. It is well-know 
that the delta of a short position in an up-and-out call option becomes arbitrarily large if, close to expiry, 
the asset price approaches the barrier. In particular this implies that A+ cannot be bounded close to expiry. 
A trader trying to hedge this position would have to buy unlimited amounts of the underlying asset. Since 
the gamma of the short position in the up-and-out call option is large and positive close to the barrier, this 



Local time and the pricing of time-dependent barrier options 



17 



delta hedge would be very profitable if the barrier were not touched. However if the barrier were broken, 
the large accumulation of the underlying asset would become a huge problem. This is why, in practice, such 
a position close to expiry would be left unhedged. Let us now proceed to the proof of Theorem l2.7l 

Proof. Let E C M+ be the finite set of discontinuities of the payoff function x ^ (l>{x)I(b-{T).b+(T)]i^)' 
which we also denote by (j) for notational convenience. We start by constructing a sequence of functions 
<pn '■ R+ K+ that satisfy the assumptions of Lemma lZjj and have the following two properties 

(1) 0„(jc) < 0„+i(jc), for allx e M+ and all n e N, and 

(2) lim„^oo (p„ (x) = ^ (x), for all x G R+ —E. 

Let p,r G E he two consecutive points in E such that p < r. In other words the function (j) is continuous 
on the interval [p,r) and has a limit at r. By Stone-Weierstrass theorem for each n G N there exists an 
element /,, e C^{[p,r]), such that the inequalities max{0(x) — ;;^,0} > /n(x) > max{^(x) — ^,0} hold for 
all X e [/?,'"). The construction implies that the sequence (/„)„gn satisfies property (1) for all x e [Pjt] and 
property (2) for all x G [p,r'). The complement R+ — E consists of a finite number of open intervals with 
the same properties as {p,r). For each point p E E we can choose a decreasing sequence of open intervals 

such that {p} = n'^^^Nl,', EHNf! = {p}, for all « e N, and A^f nA^[ = for any r e £ - {p}. Note that 
on the components of M+ — E adjacent to any p E E we have aheady constructed sequences (/„)„gn and 
(^n)nGN of functions that converge to (j) in the required way. In the complement of the neighbourhood 
NI^, we define 0„(x) to equal either /„(x) or g„{x), depending on x being larger or smaller than p. We can 
now easily extend 0„ to the interval Nn so that the resulting function is and property (1) remains true on 
any neighbourhood of p. Since {p} = n"^jA',f , property (2) is also satisfied. 

Let Z" =Z"{t,St) denote the process given in ( 12.3b that corresponds to the payoff function (p„ . It is clear 
that property (1) imphes the inequality Z"{t,S) < Z"^^{t,S) for all points {t,S) G C, where the region C is 
defined in (b) of Lemma l2.2l and all hGN. Since the set E is finite, properties (1) and (2) and the monotone 
convergence theorem imply the equality lim„^oo Z"(f ,5;) — Z{t,St), where Z, — Z{t,Si) is given by ( 12.3b . 
Since (j)„ satisfies the hypotheses of Lemma l273l expression ( 12.17b can be rewritten as 

Eo.SoMSt)]=Z"{Q,So) + 1- r A"_{t)q,{So,b^{t))dt-^ r Al{t)qt{So,b+{t))dt (2.19) 

2 Jo 2 Jo 

for all n G N and f , v e [0, T] such that t < v. For every n e N the deltas at the barriers exist by Lemma |23] 

and are given by A1{t) := limi^^ooZ'^it ,b±{t) =F e^^), where {eic)kefi is a positive monotonically decreasing 

sequence accumulating at zero. Notice that, since Z"{t,b±{t)) = for all t E [0,T], Lagrange's theorem 

implies the equalities 

for each n e N. Since the inequality Z"{t,b+{t) - £<.) < Z"+\t,b+{t) - Ek) (resp. Z"{t,b^{t) + Ek) < 
Z"+^{t,b^{t) + Ek)) holds for any n e N and all ;t e N, it follows that > AKt) > A'l+^t) (resp. < 
A"_{t) < Zi"^'(r)) for all t E [0,r]. In other words the negative sequence (zi!^(f ))„gpj (resp. positive se- 
quence {A"_ (f))«GN) is decreasing (resp. increasing) at any time t and hence converges to its infimum (resp. 
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supremum), which is not necessarily finite. We can therefore define measurable functions 4+ : [0, T] 
[-°o,0], A- : [0,T] [0,°o] in the following way A+{t) lim„^oo4!^(f), 4_(f) lim„^ocZi^(f). By ap- 
plying monotone convergence theorem to all the integrals in expression ( |2.19t we obtain formula ( 12.181 1 in 
the theorem. Furthermore, formula ( 12.18b implies that the integrals Jq A±{t)qt{So,b±{t))dt are finite. Since 
the functions A^ ,A_ do not change sign, they clearly define elements in L' ( [0, T] , ot+ (t/f ) ), L' ( [0, T] , m_ (t/f ) ) 
respectively. The system of Volterra integral equations ( 12.111 ) for the functions (4+ , 4_ ) follows in the same 
way as in the proof of Theorem |2.5l The time-dependent single-barrier case can be treated in an analogous 
way. This completes the proof. □ 

We will conclude Section |2]by considering the uniqueness of the solution of the Volterra integral equa- 
tion in ( 12.131 1 for a time-dependent single barrier in the Black-Scholes model. A much more general result 
establishing the uniqueness of the solution of the system of Volterra integral equations of the first kind given 
in (12.111 1. which requires a detailed analysis of the corresponding compact operators, will be discussed in a 
subsequent paper. 

Proposition 2.8 Assume that the payojf function (j) and the barrier b : [0, T] M+ satisfy the hypotheses of 
Theorem \2.5\ and let the asset price process {St)f^p j] follow a geometric Brownian motion. Then integral 
equation ( 12.731 ) has a unique continuous solution A : [0, T] — > R. 

Proof. Since equation (12.131 ) is linear and, by Theorem l2.5l has a continuous solution, it is enough to show 
that the only continuous solution / : [0,r] R of qu-t{b(t),b{u))f{u)du = is the obvious one, i.e. 
/ = 0. In the case of geometric Brownian motion, the integral kernel qi,^,{b{t),b{u)) is explicitly given 
by formula ( 12.81 ). By Theorem 2.1 in BTI the uniqueness of the solution of the above integral equation 
follows, if we prove that the functions k{u,t) := \Ju — tqu-tib{t),b{u)) and j^{u,t) are continuous for all 
M,f € [0, T], such that u>t, and that k{t,t) is non-zero for all f £ [0, T]. 
It follows from (12.81 ) that the function k{u,t) can be expressed as 

{B{u)-B{t)-{p-(7^/2){u-t)f\ 

2(72(M-f) ) ' 

where B{t) :— logb{t), which is clearly continuous for all m > f and has a non-zero limit, as u approaches 
f, equal to k{t,t) = This is a consequence of the Lagrange theorem (B{u) — B{t) ~ B'(^„)(m — t) 

for some G {t-u)), applied to the differentiable function B. A short calculation shows that the partial 
derivative ^ (m, f ) exists for all m > f . The regularity of the function B implies that (m, f ) has a finite limit 
at u — t and can therefore be extended to a continuous function for all t,u G [0,T], such that t < u. This 
concludes the proof of the proposition. □ 

3 Examples 

In this section we will consider some examples that illustrate the results of Theorems 12. 5l and l2. 71 We will 
look at the simplest one in Subsection |3.1| In |3.2| we solve the system of Volterra integral equations ( 12.1 It for 



k{u,t) = 
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the case of constant barriers using Laplace transforms. Subsection l3.3l brieflv discusses numerical methods 
for solving the system (12.11b in the general time-dependent case. 

3.1 Model-free barrier option price 

It is well-known that a down-and-out call option struck at K, with a barrier at the level B, has a unique 
model-independent price if B coincides with the strike K and if both the interest rates and dividend yields 
are zero. Moreover the barrier option price equals the price of a forward struck at B, as is easily seen by 
the following semi-static replication argument: since there are no interest rates or dividend yields, when 
the barrier is breached for the first time the forward contract is worth zero and can therefore be sold at no 
cost (or gain). If the barrier is not breached at all, the two payoffs clearly coincide. Since the forward has a 
model-independent price, the no-arbitrage principle implies that the barrier option price must equal S — B, 
where S is the asset price at the current time. Therefore the delta of the barrier option price is identically 
equal to one in any model. In particular the same must be true at the barrier. 

Let Ct{S,K) denote the call option price at time t in the Black-Scholes model. Equation (12.13b of Theo- 
rem l2.5l tells us that the following identity must hold 

for all t € [0, r], where the left-hand side equals the Black-Scholes formula for Ct{B,B) (the function A^(x) 
is the cumulative normal distribution) and the integrand on the right-hand side is given by ( 12.81 ) and the 
aforementioned fact A = 1. Substitution = a'^{u — t)/A and a short calculation show that this identity 
holds for all t e [0,r]. Theorem 12.51 therefore implies the following integral representation for the linear 
function S^-^ [S — B) 

S-B^C,[S,B)--j^ ^expj^ yt, 

where Cq{S,B) is the Black-Scholes formula. 



3.2 Constant barriers 

A key distinction between the constant barrier case and a time-dependent barrier case, which makes the 
former much easier to solve in a semi-analytic form, is that the kernels of the integral operators in equa- 
tions ( 12.1 lb and ( 12.13b depend only on the difference of the arguments Q{t,u) = Q{u — f ) when the barriers 
are constant. Therefore the delta along the barrier can be obtained by the following two-step procedure. 
First solve an auxiliary integral equation where the right-hand side is identically equal to one using the 
Laplace transform method, which can be applied precisely because the kernel depends on the difference of 
the arguments and the integral equation is therefore given as a convolution of two functions. In the second 
step an explicit integral representation for the delta along the barrier can be constructed using the solution of 
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the auxihary equation. In l3.2.1l we apply this method to single-barrier options in the Black-Scholes model 
(see (13.2b for the explicit formula and ll37l . Sections 8.4-1 and 8.4-4 for more details). In l3.2.2l we generalize 
this approach to the double-barrier case by finding the explicit solution of the system in (12.11b . 



3.2.1 Single-barrier options 

Let Ct{S,K) = FtN{d+) — KN{d-) denote the discounted value of the European call option in the Balck- 
Scholes model, where the forward is given by Ft :— Se^^^^'\ the drift equals 



a-.^r-d and ^± log ± 

Let a := ^^^2a^'^^ ^ positive constant and let the function q : K.+ equal q(t) :~ Let B denote 

the lower barrier {i.e.B <K) and let A :[0,T]^ M.+ be the delta of the option at level B. By Theorem l2.7l we 
need to solve the integral equation Q {B,K) = j A {u)q{u — t)du on the interval [0, T]. The substitution 
X := T — u,y:= T — t transforms the equation to 

r f{x)qiy-x)dx, ye[0,T], (3.1) 
Jo 

where f (y) := ^^Cj-yiB^K) and the unknown function / is given by /(x) :— A{T ~ x). 

As mentioned above, we first solve the auxiliary equation 1 = jQh{x)q(y'~x)dx. Recall that the Laplace 
transform of a function h is defined by ^{h){s) :— e^'"^h{x)dx for all s > such that the integral exists. 
It is obvious that ^{l){s) ^ ]■ and a short calculation yields ^{q){s) = ^s+a ' ^PPlyii^g the Laplace 
transform to both sides of the auxiliary equation we find J^{h){s) — , since the right-hand side equals 
^{h*q){s) = ^{h){s).^{q){s) by the famous property of the Laplace transform. The function {h*q){y) : = 
Jq h{x)q{y — x)dx in this formula denotes the convolution of h and q. Note that both equation (13.1b and the 
auxiliary equation can be represented in the following way: ^{y) — {f*q){y) and 1 = {h*q){y). This 
simple observation will be useful in Subsection l3.2.2l 

The task now is to compute the inverse Laplace transform which is defined as an integral along 

a path in the complex plane, of the function s i— > . Instead of using the definition of we observe 
the following elementary identities 



-V« — = and if(xi^£(vo:x))(i) 



y/s + a s^s + a s s^s + a ' 

9 r 2 

where E{x) :— [q e^*' dv is the error function. The first identity is obvious and the second follows 
from the discussion above upon noticing that the functions: E{^/ax) can be expressed as a convolution 
E{y^ax) = {y^*q){x). By applying the inverse Laplace transform to the first identity it now follows that 
the solution of the auxiliary equation is h{x) — q{x) + ^/aE{\/ax). Fubini's theorem and the auxiliary 
equation can now be used to verify that the function 



/(x):=/zWf(0) + (/i*'f")(y) 



(3.2) 
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solves integral equation ( 13.11 ). Since ^{0) = and we have the formula 

2V2 



2y/y 

the delta along the barrier A : [0, T] M+ can in the case of a down-and-out call be expressed as a convo- 
lution of two explicit functions A (f ) — {h* f ') {T — t). It is well-known that in the symmetric case when 
jJ. = 0, the down-and-out barrier option price in the Black-Scholes model is given by Co{S,K) — ^Po{S, 
where Po{S, ^) is a put option struck at ^ (see ||29l , page 454, equation (12.3)). Theorem 12.71 therefore 
yields the equation Po(5, x) = A{t)q,{S,B)dt, where q,{S,B) is given in dZSl l. 

The key observation is that the procedure described here works for an up-and-out call option in precisely 
the same way. The auxiliary equation again takes the form \ = {h-^q) (y) and hence has the same solution 
as before. Function / defined in (13.2b solves integral equation ( 13. Il l, where the function 'f' is redefined 
appropriately. In the case of an up-and-out call we have f (y) = -'^^¥.T-y,B[{ST — K)^Is^St<b}\ where the 
barrier B is larger than the strike K. In fact the same procedure works for the class of linear diffusions 
considered in Theorem 12.71 as long as one is prepared to calculate (numerically or otherwise) the inverse 
Laplace transform of the function s ^ ■^{q) {s) js. The function q: [0, T] ^ K+ would in this case depend 
on the underlying diffusion through formula (I2.7l i with x and y equal to the barrier level B. 

3.2.2 Double-barrier options 

Let B and B+ denote the lower and upper barrier respectively and let functions fi , "ft : [0,T] M be 
given by fi {y) 2(p{T ~y,B+) and ^2{y) '■= 2(p{T —y,B^), where (p represents the discounted value of 
the European payoff (see Theorem l2.5l l. By introducing the change of variable x :— T — u as in the previous 
subsection and denoting f\ (x) := A+{T — x), f2{x) := A^{T —x), we can express equation ( 12.1 11 1 using the 
linear operator .yT : ( [0, T] ) x l' ( [0, T] ) ^ ( [0, T] ) x l' ( [0, T] ) in the following way 

where ^ f/^V,) f/o'en + /^'en(y-.)/2 

Vf2y [fij [fij \I^Q2i{y-x)Mx)dx + syQ22{y~x)f2{x)dxJ 

The functions Qij : [0, T] — > M, y e { 1 , 2}, are the coordinates of the matrix in ( |2.12| i and can be expressed 
as functions of one variable precisely because the barriers are constant in time. 

Recall that convolution can be used to make the Banach space L'([0,r]) into a commutative Banach 
algebra, since the function {u*v){y) — Jq u{y~x)v{x)dx is an element of L' ([0, T]) for any u,v E {[0,T]). 
Using this multiplicative structure and the definition in (I3.3l l we can express the linear operator J(f in the 
following way 

^fi\ fQnQn\Jfi' 
,G2i 622/ * V/2 

The integral equation in ( 13.3b can now be solved in two steps. The first step consists of finding the 
functions hij : [Q,T] R, i,j e {1,2}, which satisfy the identity 

621 622 / \h21 h22 / 
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Note that the product of any pair of coordinate functions in this expression is given by their convolution. 
Note also that the solution of this auxiliary equation depends solely on the barrier levels B+,B and is 
independent of the payoff of the option we are trying to price. By applying the Laplace transform ^ to 
each coordinate of this equation, we obtain a linear system for the functions ^{hij), where multiplication 
is defined using a point-wise product rule: 

i 0\ ^ (^{Qn){s) ^{Qi2){s)\ (^{hn){s) ^{hn){s)\ 
\^iQ2i){s) ^(Qiim) [^{h2iKs) ^{h22){s)l ■ 

Assuming that the determinant {.^{Qiij-^iQii) — -^(2i2)-^(G2i))(*) is non-zero for all s > 0, we can 
explicitly solve this system of equations. In order to obtain the functions hjj : [0, T] M, /, /' G {1,2}, we 
need to perform Laplace inversion on each of the four coordinates of the solution of the linear system. 

Once the auxiliary equation in ( 13.4b has been solved, we can express the solution of the original integral 
equation in i3.3i in the following way 

Mx)\ ^ fhnix) hn{x)\ /fi(0)\ ^ fhn hn\ ^ H\ 

/2 W ; \h2l (x) h22ix) J Vf2(0) ; \h21 h22 J ^ \%) 

Since in our case we have fi (0) = ¥^(0) = 0, the deltas at the upper and lower barriers are given by the 
formulae A+{t) = (/in * W[){T - r) + [hn * %){T - 1) and A^{t) = {h2i * fi')(r - f ) + * I'i^T - 1) 
respectively. Representation ( 12.181 ) of the double-barrier option price in Theorem l2.7l can now be applied. 

3.3 Time-dependent barrier options 

In case of general time-dependent barriers not much can be said analytically about the structure of the 
solutions of the system of integral equations in ( 12.111 ). However the trapezoidal product integration method, 
described in BTI . can be applied directly to the single-barrier problem. The substitutions y := T — t and 
X :— T — u, used in Subsection 13.21 transform equation (12.13b into a generalised Abel equation with a 
weakly singular kernel 

k{y,x):=T^^^qy-Ab±iT-y),b±{T-x)), 

where the function b± : [0, T] is either a lower or an upper barrier and the function q is given in ( 12. 7b . 

Using the notation ^{y) := (p{T —y,b±{T—y)) for the discounted value of the European payoff (see The- 
orem |23] for the precise definition of (p) and f(x) := A±{T — x) for the unknown function in our integral 
equation, we can rewrite ( 12.13b as follows: 

W{y)^ f ^^^f{x)dx. (3.5) 

It follows from the representation of the density of a linear diffusion given in fT5l|, Section 1.2, equa- 
tion (2.8), that k{y^y) := lim;^/^ k{y,x) exists and is non-zero since we are assuming that the barrier function 
is differentiable. This statement is clear for geometric Brownian motion, when the function q is given 
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by (I2.8I 1, and essentially the same proof can be used for a general linear diffusion once we apply the rep- 
resentation given in ifTsll . The observation that < k{y,y) < 00, for all y £ [0,T], is of utmost importance 
because it makes the lower-triangular system of linear equations in BTl non-singular 

The main theorem of [41] says that the solution of the lower-triangular linear system, given by equa- 
tion (4.1) on page 179 of the same paper, converges to the solution of the integral equation i3.5i at the 
order of 0{h^), where h is the distance between the consecutive points in the discretization of [0, T]. This 
convergence result assumes some regularity properties of the solution /, such as continuity on the entire 
interval [0,r], which are in general not satisfied in our context. The situation is improved if we work in 
the domain of Theorem l2.5l In other words when faced with a discontinuous payoff function ^, we can 
first approximate it by a smooth function 0„, as in the proof of Theorem l2.7l and then solve the linear sys- 
tem from [41 1 which corresponds to the derivative that delivers This procedure introduces an additional 
numerical error since we are pricing the "wrong" derivative, but improves the convergence speed of the 
solution of the lower-triangular linear system from PTI . It follows from the construction of 0„ in the proof 
of Theorem 122] that the price of the barrier option with the payoff (/)„ converges uniformly in (f ,5,) to the 
price of the same barrier option with the payoff ^. The stability of the proposed numerical algorithm will 
be the subject of future research. 

The double-barrier case can be dealt with similarly upon noticing that the functions qu-t{b+{t),b^{u)) 
and qii^t{b-{t),b+{u)), which appear "off the diagonal" in the kernel of the system of Volterra equations 
given in (12.121 1. are smooth and bounded for all f , m e [0, T] such that f < m. In other words we can extend the 
n-dimensional lower-triangular system from BTl . used to solve integral equation ( 13.5b . to a 2«-dimensional 
lineal" system by representing the integrals against functions qu-t{b+{t)^b^{u)) and qu-t{b-{t),b+{u)) us- 
ing the standard trapezoidal method (which can be expressed as matrix vector multiplication). By express- 
ing the solution vector in the following way (zi+(fi),Zi_(fi), . . . ,Zi+(f„),Zi_(f„))^, where (f,),= i „ is an 
increasing sequence such that f 1 = and t„ = T, the 2n-dimensional linear system we need to solve becomes 
lower-triangular because of the identities lim„\j qu-t{b+{t) ,b-{u)) = limii\t qu-t{b-{t),b+{u)) = 0. There 
are a number of algorithms designed to solve this kind of linear system very quickly and accurately (see for 
example [28 J). Their implementations usually rely on numerical libraries like BLAS and LAPACK for the 
calculations. These numerical libraries are highly optimised and can be called directly from C++ (see ^ 
for more information on LAPACK, also available at http://www.netUb.org/lapack/lug/). However the ques- 
tion of implementation and the optimal choice of algorithm for solving our lower-triangular linear systems 
requires further numerical investigation. 

4 Conclusion 

In this paper we have obtained an integral representation of the difference between the time-dependent 
double-barrier option price and the price of a European option with the same payoff. Theorems I2.5l and l2.7l 
give the precise formulae in terms of the double-barrier option deltas : [0,r] ^ R x K at the 

barriers (see (a) in Lemma |23T l for the precise definitions of these functions), which solve the system of 
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Volterra integral equations of the first kind in ( 12.1 lb . It follows by construction that the system of integral 
equations in (12.111 1 has a solution. The most natural question is the one of uniqueness of solution, which is 
equivalent to the question of whether zero is an eigenvalue of the compact linear integral operator with a 
weakly singular kernel given in (I2.12l i. In general, compact operators exhibit both kinds of behaviour and 
the standard technique of transforming an integral equation of the first kind Kf = 'f' to an integral equation 
of the second kind / — L f = f , where K and L are related integral operators (see [37J, Section 8.3), does 
not apply in our setting because of the weak singularity of our kernel. The transformation of the problem is 
desirable because we can use the Fredholm alternative to analyse the kernel of the operator I L, where L is 
compact and / is the identity operator. In spite of this difficulty this general approach can be made to work 
in the case of equation ( 12.1 lb by careful inspection of the construction of the functions in the integral kernel 
and the uniqueness of the solution can be proved, at least for payoff functions that satisfy the regularity 
conditions of Lemma l23] The proof will be given in a subsequent publication. 

If we assume that the Volterra integral equation of the first kind in (12.1 lb has a unique solution, then 
representation ( 12.10b for the time-dependent double-barrier option price from Theorem l2.5l implies that the 
value of this path-dependent derivative depends only on the one-dimensional distributions in the risk-neutral 
measure of the underlying process 5,. It is well-known that having the vanilla option prices for all strikes 
and all maturities is equivalent to having the one-dimensional distributions of St- Theorem 12.51 therefore 
provides an explicit link between the vanilla option prices and the barrier option prices for all reasonably 
smooth barriers and a wide class of local volatility models. The fact that vanilla option prices determine 
uniquely the barrier option prices in the world of local volatility models has been known since the seminal 
work of Bruno Dupire 1 1 2 1 where, under certain regularity conditions, a PDE for the local volatility function 
X I— > a{x) is derived from the vanilla option prices. This in turn determines the risk-neutral dynamics of St 
and therefore the prices of all path-dependent derivatives but does not yield an explicit relationship. 

In this paper we have discussed a barrier pricing problem without rebates. In general a barrier can pay 
a contract defined rebate F+ (f ) (resp. F-{t)) if at time t E (0, T) the asset price process Si equals the barrier 
level b+{t) (resp. (f )) for the first time since inception. It is not difficult to see that under some additional 
technical assumptions on the functions F± : [0, T] R, the change-of-variable formula from (35] can be 
applied and a similar technique to the one used to prove Theorem l2 . 5 1 vields an integral representation of the 
time-dependent double-barrier option price with rebate. In fact the final formula is very similar to the one 
in ( 12.101 ). with (p{Q,So) replaced by the sum of the expectation of the European payoff and certain integrals 
over the time-interval [0, T] of the rebate functions F±. 

We have seen that purely probabilistic concepts such as local time and the generalized Ito formula 
proved by Peskir in |35 1 can be used to obtain a new structure in the barrier option pricing problem, which 
can then be applied to the pricing and hedging of double-barrier options in local volatility models. This 
structure consists of two deterministic functions A± : [0, T] R, which represent the deltas at the two 
barriers. It is intuitively clear that the same structure exists in stochastic volatility models. Since Peskir's 
formula has been generalized to higher dimensions in 1361 for all (possibly discontinuous) semimartingales, 
a generalisation of the approach presented here might be feasible. The multidimensional change-of-variable 
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formula in ll36l is both surprising and satisfactory, not just because it applies to all semimartingales but 
because, under natural deterministic conditions on the value function, the resulting formula is a direct 
extension of the one-dimensional formula in |351. However a direct application of the formula in ll36l is not 
possible in our case, because it requires the existence of a regular extension of the value function across the 
boundaries of its natural domain, a question that requires some investigation. The issue of which quantity 
one could represent in terms of deltas at the barriers in the higher dimensional case (the value function itself 
cannot be represented) provides in our view an additional interesting problem for future research. 



Appendix 

A A change-of- variable formula with local time on curves 

In this section we establish a mild generalisation of the change-of-variable formula given in Theorem 3.1 
of |[35l . In fact Theorem lA.ll is implicitly proved in ||35l . Since Theorem lA.ll is central to our analysis, for 
complicity we give a proof based on a direct application of Theorem 3.1 and Remark 2.5 in ll35l . 

Let X := (^f)fG[o.7-] be an Ito diffusion that solves the following stochastic differential equation dX, = 
M{Xt)dt + E{Xt)dWt, where M{x) := jxx and E{x) -.^ xa{x), and let b± : [0,T] ^ M be two continuous 
functions offinite variation satisfying Z7_(f) < Z7+(r) for all f e [0,r]. As before we set C := {(r,x) G [Q,T)x 
M; b-{t) <x<b+{t)},B+ := {{t,x) e [0,T) xR; x > b+{t)} and B := {(f,x) e [0,T) xR; x < b-{t)}. 
Let F : [0,r] X Mbe a continuous function which is C''^ on the open subset B_ UCUB+ of [0, T] x R. 
Given a function g : [0, T] ^ K of bounded variation let V (g) (f ) denote the total variation of g on [0, t] for 
any t <T . 

Theorem A.l Let T,{x) > Q for all x € (0,oo) such that {t,x) £ C and assume that Ft + MF^ + ^F^x is 
locally bounded on B_ U C U B+, the limit Fx{s, b± (s)±) := lim^^^o Fx{s, b± (s) ± e) is uniform in s S [0, f] 
and that swpQ^g^gV{F{-,b±{-)±£)){t) <°°forsome5 > Q and any combination of sings + and — . Then 
the following change-of-variable formula holds: 




For the definition of the local time of X at the curve b, ll^iX), see page|8] 
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Proof. We start by defining continuous functions : [0, T] x M ^ M which satisfy the hypothesis of 
Theorem 3.1 in 1351 and then apply the theorem to obtain the formula above. Since functions b± are con- 
tinuous, the images b±{[0,T]) are disjoint compact subsets in with strictly positive distance. Hence 
there exist e > such that 4e < b+{t) — b^{t) for all f e [0, T]. It is clear that there exist smooth functions 
c± : [0,r] ~> M that satisfy b^{t) < c_(f) < b-{t) + e and b+{t)-e < c+(f) < b+{t) for all f e [0,r]. In 
particular it follows that c+(f) — c_(f) > 2e. 

We now define continuous functions and F , which are C''^ everywhere in [0, T] x R except along 
the curves b+ and b- respectivly, by the formulae 

jFit,x) ifx>c_(f), jF{t,x) if x<c+{t), 

F [t,x) :— < F (t,x) :— < 

[ f^{t,x) ifx< c^{t), y .nt,x) ifx > c+(f). 

The function /+ (resp. /^) is a C''^ extension of F across the smooth boundary c_ (resp. c+), which exists 
because F is C''^ by assumption on the domain C. Note also that the functions are non-unique. Given 
these definitions of F*, the only discontinuities of the derivatives are the ones inherited from the original 
function F along the curves b± respectively. Since F satisfies the conditions in Theorem lA.il the functions 
also satisfy the assumptions of Theorem 3.1 in ||35]| . which therefore implies the following formulae for 
any fixed time t £ [0,T]: 

F^{t,X,) = F±(0,Xo) + ^'(F± +MF± + ^F,^^){s,X,)I{x^^,^^,^}ds + J\EF,^)is,^^^^^^ 

where the signs ± are simultaneously equal to either + or — . 

Let {^t)te[a.T] denote the filtration of the Brownian motion {Wt)te[o.T] th^t satisfies the usual conditions. 
Since the processes on both sides of the equality in ( lA.lb have continuous paths we can assume that they are 
indistinguishable and therefore substitute fixed time t e [0, T] with any stopping time relative to )te[Q.T] ■ 

We now define an increasing sequence of stopping times in the following way: p\ t /\i\\f{s;Xs = 
c+{s)}, p2 f Ainf{i > p^X, = c_(i)} and p2„+i f Ainf{.s > P2„;X, = c+(i)}, p2„+2 f A inf{i > 
P2n+i;^i = C-(j')} (we are assuming wlog that Xq < c+(0)). Note that for any s in [p2ii+i,P2n+2] (resp. in 
[p2„,P2,7+i]) the value of the random variable Xs is strictly above b-{s) (resp. below b+{s)) and therefore 
F{s,Xs) equals F^{s,Xs) (resp. F^{s,Xs)). We also have lim„^ocP„ = t a.s. and, for almost all paths of X, 
p„=t for some n g N (this follows from the inequality c+ {s) — c_ {s) > 2e, for all s G [0,T], and the fact 
that the expectation of the upcrossing number of our semimartingale is finite, cf. Il25l . Theorem 1.3.8 (iii)). 
For a fixed f e [0, T] we have a telescoping representation 

F{t,X,)-F{0,Xo) = £ {Fip2n+2,Xp,,^^2)-F{p2n+uXp^„^^)+F{p2„+uXp^„^,)-F{p2n,XpJ) 
n=0 

= {F^ {P2n+2,Xp2„^2) ^ F^ iP2n+uXp2„+i)) + (■f'" (P2«+l ,^P2„+i ) ^ (P2,n^P2„ ) ) i 
n=0 n=0 
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where po := 0. We are allowed to reshuffle the summands in this path-wise identity since for almost all 
paths the sums consist of finitely many summands. The theorem now follows by applying formula (lA.U 
between the stopping times for (^,) to the summands in the last expression and collecting the terms. □ 

B Analyticity properties of time-dependent barrier option prices 

Let (p : M+ M.+ be a payoff function that is continuous on a complement of a finite set where it is 
right-continuous with left limits. In particular (p is continuous at zero. An important example is ^{x) = 
(jc — g^)(x) for some constants B- <K < B+ wdx e K+. The diffusion X = {^i)t£[o.TY specified by 
the time-homogeneous SDE c/X; —M{Xt)dt+'L{Xt)dWt with linear drift M(^) :— /ijc and a locally Lipschitz 
diffusion coefficient Z (x) :^xa{x) > 0, for x G (0,°°), is as described in the beginning of Section|2] Using 
the notation from appendix lAl we define a family of stopping times T,, for any f g [0,r], by T, := inf{v €E 
[0,r — f]; e M— (^7_(f + v),^7+(f + v))}, where the boundary functions /j-i- : [0,T] M are continuous 
and twice-differentiable in the interval (0,r). We also consider the case where either b+ or are not 
present to capture the single-barrier case. If there is no upper barrier we assume in addition that (j){XT) G 
L'(i2,Q). The discounted barrier price process is a martingale given by := E [0(XTo)/{To=7'}l-^r] where 
the filtration {•^t)te[o.T] is as described in the beginning of Section|2] The identity /{^p^T-} = 
the facts {tq > f } G Tq = f + T, on the set {Tq > ?} and the Markov property of X imply the following 
path-wise representation for the barrier price Vt ~ /{T(,>r}Z(f where the function Z : [0, T] x R ^ M is 
given by 

Z(f,x) := E,., [0{Xr+r,)I{r+r,=T}] ■ (B-D 

Here the process X starts at time t with value X, = x. 

It is often stated that the barrier option price satisfies a certain PDE with absorbing boundary conditions. 
Such statements are in fact referring to the analyticity properties of the function Z, which we make precise 
and prove in Theorem lB.il 

Theorem B.l Let ^g{t,x) := (^gt +Mgx+ ^gxx^ (Itx) be the infinitesimal generator of the diffusion 
^ — (^?)rG[0,r]> where Yt :— {t,Xt) and the process X is as described above. Let the set C be as defined 
in appendix^and assume that C ( the closure is taken in the space M x [0, T]) is contained in (0, oo) x [0, T]. 
Then, under the above hypothesis on the payoff (p and barriers b±, the function Z given by HB.l\l is contin- 
uous on the set C — (M x {T}) and solves the following parabolic boundary value problem: 

^Z{t,x) for (t,x) gC, 
Z(r,x) = 0(x) for X G {b-{T)M{T)), 
Z{t,b±{t))=QfortG[QJ]. 

In the single-barrier case, the local behaviour and the terminal condition satisfied by the function Z remain 
the same, but the boundary conditions change as follows: for an up-and-out option the boundary conditions 
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are Z{t,b^{t)) = and Z{t,0) = ^(0), t G [0, T], and in a down-and-out case we have Z{t,b- (t)) ~ 0/or 
allte [0,T]. 

Proof. Assume first that both barriers are present. If in addition we assume that (p : [b^{T),b+{T)] R 
is continuous and satisfies (p{b-{T)) — ^{b+{T)) — 0, then by Theorems 6.3.6 and 6.5.2 in lfT6l . function 
Z defined by dB.lb is the solution of the parabolic PDE and satisfies the required boundary conditions. In 
particular Z is continuous on C. 

Assume now that the payoff is discontinuous. By assumption ^ only has a finite number of bounded 
jumps . Hence we can express it as ^ = lim„^oo 0„ , where the functions 0„ are continuous with (j)„{b^{T)) = 
^n{b+{T)) = for all « G N and the convergence is uniform on the complement of any neighbourhood of the 
discontinuities of 0. In fact we can choose the functions 0„ so that there exists a decreasing sequence of open 
sets A^„ C M, such that the intersection jA^„ equals the set of discontinuities of 0, and {x) = (/)„ [x) on the 
complement of N„ for all n e N. In the obvious notation we get |Z(f ,x) — Z„{t,x)\ < AKt,x[lN„ {Xt)], where 
A is some constant independent of n which exists since ^ is a bounded function. As usual I^,, denotes the 
indicator function of the set Nn ■ Since the random variable Xj has a density in the set (0, °o) which is smooth 
in the parameter (f ,x) (see the discussion preceding Theorem l2.5l in Section|2|l, the right-hand side of the last 
inequality goes to zero uniformly on some neighbourhood of the point {t,x). This implies that Z is a limit 
of a uniformly convergent sequence of continuous functions and is therefore continuous on the complement 
of the finite set of discontinuities of the payoff (j). Note that Z{t,b±{t)) = lim„^oo Z„{t,b±{t)) = for all 
t e [0, T) and that Z{T,x) = 0(x) by definition. 

We now need to prove that Z is in C^'^(C) and that it satisfies the PDE ^Z = 0, where is the 
infinitesimal generator of the diffusion Y. These are local properties of the function Z and it is therefore 
enough to show that they hold on any bounded neighbourhood U C C of an arbitrary point {t,x) E C. 
We can assume without loss of generality that the boundary dU is smooth. Then the parabolic boundary 
value problem for g : U ^ R, given by Jifg = in U and g\gij — Z\gij, has a unique solution (see fT6l . 
Theorem 6.3.6). Let be the first time the process Y, which started at (f ,x) e U, hits dU. From Dynkin's 
formula (see ll33l . Theorem 7.4.1) and the fact that g is the solution of the above Dirichlet problem we find 
g{t,x) = E( ;t[Z(Fi+T(/)]- Since Z satisfies the mean-value property (see |f33l . page 121, formula (7.2.9)), it 
follows that g{t,x) = Z{t,x) for all {t,x) S U. This proves the theorem in the double-barrier case. 

If we only have a lower barrier, we can express the function Z as a limit of double-barrier option prices 
where the "artificial" upper barrier tends to infinity. Using a similar argument as above, and the fact that the 
maximum of the process X is finite Q-almost surely in the time interval [0, T] (by Lemma lTTl the process X 
does not explode to infinity in finite time) it is not hard to see that the convergence is locally uniform, which 
in turn implies that the function Z is continuous on the complement (in C) of the discontinuities of (j). Once 
we have established continuity, the same "local" argument as in the paragraph above proves the theorem in 
the case where there is no upper barrier. 

In the up-and-out case, we introduce a constant lower barrier at some small level e with the boundary 
condition Z(f , e) — 0(e) for all t e [0, T]. Since the function (j) is continuous at zero, the argument similar 



Local time and the pricing of time-dependent barrier options 



29 



to the one above yields continuity of the solution of the parabolic problem obtained in the limit as £ ^ 0. 
Once we have continuity of the solution, the "local" behaviour follows as in the preceding two cases. □ 



C Proofs of Lemmas|2ll|23]and|22] 



We start by proving Lemma 1271] 

Proof. By Feller's test for explosions (see Theorem 5.5.29 in ||25]| ) it is well know that the statement of the 
lemma holds if the following iterated integral diverges 

iz X 

for some z,xq e (0,°°). This is clearly true if the limit \\m.x^„B{x) is finite. If this is not the case, a simple 
application of L'Hopital's rule implies that the functions ^ exp(B(x)) /^°°exp(— is asymptotically 
equal to the function x jca(x)^. This proves the lemma. □ 



Next is the proof of Lemma 1231 

Proof. The lemma is a consequence of Schauder's boundary estimates for the solutions of the initial 
parabohc partial differential equations proved in [14]. Let us first consider the double-barrier case. Recall 
that C denotes the domain of the solution of the PDE from (b) of Lemma IZ21 Denote by F : dC — ({0} x 
(0) ,b+{0))) M+ a continuous function which maps the curves b±{[0,T]) to zero and coincides with 
the payoff on the interval [b^{T),b+{T)]. By Theorem 3.3.7 on page 65 in ifTSl the partial derivatives 
Zs , Zss , Zt of the solution Z of the PDE in (b) of Lemma l272l will be Holder continuous of order a G (0 , 1 ) on 
C, if we can find an extension f : C — > K+ of the function F, whose partial derivatives "f^, f^s, "f? are Holder 
continuous of order a on the domain C. Theorem 3.3.7 in |,15J applies in our case because the volatility 
function x i— > xa{x) is uniformly elliptic on the domain C since it is strictly positive on the compact set C. 
Before constructing an extension W explicitly, let us show that Holder continuity of the partial derivatives 
of Z implies the lemma. 

Pick a sequence (e„)„GN of positive real numbers which converges to zero. Since the second derivative 
Zss is Holder continuous on a bounded domain C, its modulus must be bounded by some constant c (i.e. 
|Z55(/,;ic)| < c for all points {t,x) e C). Therefore, by Lagrange's theorem, we have \Zs{t,b+{t) — e„) — 
Zs{t,b+{t) — £k)\ < c\e„ — Ekl for all n,k e N and all t G [0,T]. Since the right-hand side of this inequality 
is independent of time f, the sequence of functions (f k-> Zs[t,b+{t) — e„))„GN is uniformly Cauchy on the 
interval [Q.T] and therefore converges uniformly to the continuous limit A^. The same argument can be 
used for the lower barrier. This implies part (a) of the lemma. 

For part (b) let us choose a real number 5 > such that the point {t,b+{t) — 25) lies in the domain C for 
all t G [0, T]. Since the barrier b+ is uniformly continuous on the interval [0, T], there exists 5q>0 with the 
following property: if |f — i| < 5o, then \b+{t) — b^{s)\ < 5 for all s,t S [0, T]. Choose any e e (0, 5) and 
assume that f,- G [0, T] satisfy = fo < f i < ■ ■ ■ < t„ ~ T and max{f,- — f,_ i ; / = !,...,«} < Sq. Note that this 
implies that, if b+ {ti) > b+ {ti^ i ), for any / e { 1 , . . . , n — 1 }, the point , b+ {ti^ i ) — e) lies in the domain C. 
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Similarly if < /7+(f,_i) we find that the point {ti^\,b+{ti) — e) is in C. Using these observations we 

obtain 

\Z{ti,b+ (fi) - e) - Z(f;_i , b+ (f;_ 1 ) - e) I < \Z{ti,b+ {ti) - e) - Z(fi_ i , ^7+ (f,-) - e) | + 

|Z(f,-_ 1 , /7+ (f,-) - e) - Z(f,-_ 1 , b+ (f,_ 1 ) - e) I 

< A (f - f 1 ) + D I /7+ (f ) - ^+ (f/- 1 ) I 

< (A+D£)(f,--f,_i), 

where the constants A, D,/? are upper bounds on the absolute values of the derivatives ^jZ^,/?'^ respectively. 
In this inequality we assumed that < /7+(f,_i). In case b+{tj) > i) a similar bound with the 

same constants, which is also independent of e, can be obtained. This inequality implies that the family of 
functions (f t-^ Z{t,b+{t) — £))ee{Q.S) has a uniformly bounded total variation. Lower barrier can be dealt 
with in an analogous way. This proves part (b) of our lemma. 

We are now left with the task of showing that the payoff (p can be extended to a function ^ .C^ M+ with 
Holder continuous derivatives Wt,^s,^ss- We start by defining a global diffeomorphism j3 : [0, T] x R — > 
[0, T] X M, which straightens the barriers of the region C, given by 

/3(r,5):-(r,B(r,5)), where Bit,S) :^^±^j)-^-^{S-b4t))+b4T). 

b+{t)-b^[t) 

Note that P{t,b±{t)) = (t,b±{T)) for all f e [Q,T] and I5{T,S) = {T,S) for all 5 G R. We can therefore 
define ^{t,S) := (p{B{t,S)) for any point {t,S) EC. A simple calculation shows the following: 

b+{T)-b^{T) 



b+{t)-b^{t) ' 
b+{T)-b^{T) 



^sit,S) = f (B(f,5)) 

The desired properties of the function f follow directly from the assumptions in the lemma on the payoff 
and the boundary functions b± . 

Our final task is to prove the lemma in the case where there is only one barrier. Theorem 3.3.7 on page 65 
in ifTSl can only be applied if the domain C is bounded. Assume we only have, say, a lower barrier f t-^ b{t). 
Then by Theorem lB. li the discounted time-dependent single-barrier option price Z{t,S) still solves the PDE 

from (b) of Lemma 12.21 We can now introduce artificially a constant upper barrier at some large value 

S~ (J^ (s) 

B and formulate a parabolic initial-boundary value problem Ut{t,S) + jiSUs{t,S) H 2^^ssit iS) — 

on the bounded domain C -.^ {{t,S) : t e [0,T],S e [b{t),B]} with the payoff function (p : [b{T),B] ^ R 
and boundary conditions U{t,b{t)) ~Q,U{t,B) = Z{t,B) for all t G [0, T]. Like in the double-barrier case, 
because the domain C' is bounded, our assumption on the volatility function ff implies that the differential 
operator is uniformly elliptic. By Theorem 4 in [ L4J such a problem has a unique solution and therefore 
U{t,S) — Z{t,S) for all {t,S) e C. Furthermore the same argument as above imphes that functions 1 1-^ 
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Us{t,b{t) + e) converge uniformly to a continuous function 1 A{t) defined on [0,T] and that the total 
variation of 1 U {t ,b{t) + e) is bounded uniformly for all small positive e. This proves the lemma in the 
case of a single-barrier option with a lower barrier. The single upper barrier case can be dealt with similarly. 

□ 

Finally we demonstrate Lemma l2!6l 

Proof. Since the set K is compact, it follows from definition ( 12.11 ) that it is enough to prove the lemma for 
the transition density p{u — t;x,y) of the process S. As mentioned on page [10] (see also 1201 . page 149) 
the function (m,x) — > p{u — t;x,y) solves the parabolic PDE problem on the bounded domain {t,T] x K. 
Let the function (m,x) — > v{u — t;x,y) be the solution of the same PDE satisfying the boundary conditions 
v(m — t;x,y) = —p{u ~ t\x,y) for all m G [f , T] and x e dK (the symbol dK denotes the two boundary points 
of the interval K) and the initial condition v{Q\x,y) = for all x e — dK. Such a solution exists and is 
bounded because the PDE is uniformly parabolic on the domain (f,r] x K and the boundary conditions 
are continuous and bounded. Furthermore, by Section 5 in ifTSIl . there exists a non-negative fundamental 
solution / for our parabolic PDE that satisfies the inequahty in Lemma 12.61 By the maximum principle 
(see ifTSl . Theorem 2.1 on page 34), the function / dominates the solution v + p of the PDE on the entire 
domain (f , T] x K. Since the v is bounded the lemma follows. □ 
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